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ABSTRACT 

We use non-linear scaling relations (NSRs) to investigate the effects arising from the 
existence of negative correlations on the evolution of gravitational clustering in an 
expanding universe. It turns out that such anti-correlated regions have important 
dynamical effects on all scales. In particular, the mere existence of negative values for 
the linear two-point correlation function fj. over some range of scales starting from 
I = L , implies that the non- linear correlation function is bounded from above at all 
scales x < L a . This also results in the relation £ cx x~ 3 , at these scales, at late times, 
independent of the original form of the correlation function. Current observations do 
not rule out the existence of negative £ for 200ft- -1 Mpc £ 1000/i _1 Mpc; the 
present work may thus have relevance for the real Universe. The only assumption 
made in the analysis is the existence of NSR; the results are independent of the form 
of the NSR as well as of the stable clustering hypothesis. 
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1 INTRODUCTION : THE IMPORTANCE OF 
BEING NEGATIVE 

The problem of gravitational clustering in an expanding Uni- 
verse is one of the important open issues in cosmology. In the 
standard picture, non-linear structures were formed through 
the mechanism of gravitational instability, via the amplifica- 
tion of small initial density perturbations. At early epochs, 
when the fluctuations are in the linear regime, perturba- 
tion theory can be used to study the evolution. Analytical 
understanding of the quasi- and non-linear regimes has un- 
fortunately proved a more intractable problem and much of 
the activity in this area has centred on numerical simula- 
tions. It is, however, important to understand the numeri- 
cal results; a number of attempts have hence been made to 
address the subject by semi-analytic means. While this pro- 
gram is still in its infancy and we are far from providing a 
"first principles" paradigm, such attempts have nevertheless 
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been successful in identifying at least some of the key issues 
in the study of gravitational clustering. 

A key problem in this area involves understanding how 
power gets transferred between different scales during the 
evolution of clustering. In the simplest context, one would 
like to understand how an initial power spectrum, peaked 
aroun d a particular scale, evolves w ith time. Numerical re- 
sults (Bagla & Padmanabhan 1997) and a few analytic in- 
sights (especially regarding inverse cascade; see Padmanab- 
han 2000) show that power injected at a scale L cascades 
both downwards (to x < L) and upwards (to x > L). To 
the lowest order of approximation, the Fourier transform 
of the gravitational p otential is described by the equation 
dPadmanabhan 2000|): 
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which governs the growth of 0k due to non-linear mode cou- 
pling. When the non-linear coupling term on the right hand 
side is small, we get the standard result that 0k is inde- 
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pendent of time in the linear regime. If, however, the ini- 
tial power spectrum is sharply peaked at some scale, the 
issue becomes more subtle. Linear theory remains valid at 
all scales only as long as 0k is small at all scales. As soon 
as some scale goes non-linear, it is possible for the integral 
on the right hand side of equation (Q) to pick up significant 
contributions to drive the evolution of 0k at other scales. For 
example, this will lead to a fc 4 tail in the p ower spectrum 
of the density contrast ( Padmanabhan 2000| ) , purely due to 
non-linear coupling. This is a clear case of inverse cascade 
of power in gravitational clustering. 

Power transfer also occurs to smaller scales but is some- 
what more difficult to analyse from first principles. Simula- 
tions ( Bagla fc Padmanabhan 1997) as well as theoretical 
arguments (Padmanabhan & Engineer 1998) suggest that 
the evolution to smaller scales will lead to a universal form 
of the power spectrum, reminiscent of the Kolmogorov spec- 
trum in fluid turbulence, if the initial spectrum was sharply 
peaked at some scale. 

The above results emphasise the simplicity and tractability 
of a question involving the evolution of a power spectrum 
which was originally peaked around some scale. This issue is 
important in understanding gravitational clustering, in spite 
of the fact that currently fashionable models for structure 
formation have broad band power spectra, with P(k) oc fc for 
small fc. To study the evolution of such power spectra ana- 
lytically, we need to use some technique which is sufficiently 
general and yet tractable. One such approach is based on 
non- linear scaling relations (NSRs, hereafter), which provide 
a mapping between th e linear and non- linear mean two-point 
correlation functions (Hamilton et al. 1991; Nityananda & 
Padma: labhan 1994; Mo at al. 1995; Padmanabhan 1996a). 
The above mapping appears to be validated by numerical 
studies, which have also shown it to be fairly model inde- 
pendent. Given an initial mean correlation function, £l, in 
the linear regime, the NSR yields the evolved mean correla- 
tion function, £, (and thus, the power spectrum) at any later 
time. It also illustrates how power is transferred between dif- 
ferent scales during the evolution of clustering. Finally, the 
limiting forms of the NSR can be derived using the standard 
paradigms of scale-invariant radial collapse and stable clus- 
t ering in the quasi-line ar and non-linear regimes respectively 
( Padmanabhan 1996a| ); it has also been shown that scaling 
laws are a generic f eature of clustering in an arbitra ry num- 
ber of dimensions (Padmanabhan & Kanekar 200C). 



Non-linear scaling relations have thus proved a power- 
ful tool for understanding gravitational clustering; however. 



they h ave one major defect. The standard NSR (Hamilton 
et al. 1)91), by its very form, is restricted in its applicabil- 



ity to positive semi-definite correlation functions; negative 
values of £ are not included in its gamut. However, any suffi- 
ciently localised power spectrum (for example, a Dirac delta 
function), will produce negative £ at some scales; in fact, 
as will be seen later, negative £ can also arise from broad 
band power spectra. It is thus interesting, from a purely 



mathematical perspective, to investigate the effects of anti- 
correlation on the formation of structure. 

The standard NSR shows that power is transferred from 
large scales to small scales during the evolution. It turns out, 
however, that the very existence of negative £, combined 
with the equation of conservation of pairs, implies that the 
flow of power will reverse direction at the scale at which £ 
crosses zero. Further, it will be shown that the transfer of 
power to small scales saturates at late times; this results in 
an upper bound on £ at all scales smaller than the one at 
which £ is initially negative. Finally, the correlation function 
at these scales satisfies the relation £ oc x~ 3 , at late times, 
independent of its original form. 

None of the above features arise for initial power spec- 
tra with positive definite correlation functions. It should be 
stressed, however, that these results are valid even if ^ is only 
negative over a very small range of scales or at arbitrarily 
large scales; the results also do not depend on the extent to 
which f goes negative but arise entirely due to the existence 
of negative £. They are also quite independent of the sta- 
ble clustering ansatz. Current observations do not preclude 
negative values of the correlation function at some scales; in 
fact, we will show in section ^ that one can indeed modify 
the power spectrum of the universe, so as to have negative 
correlation functions over certain ranges of scales, without 
violating any observational constraints. 

The plan of the paper is as follows : the mathemat- 
ical framework is set up in section ^; we also discuss the 
usual NSR for the positive £ regime here. Next, NSRs for 
other statistical indicators of clustering, such as the power 
spectrum, are considered in section ^. We move on to re- 
view observational constraints on the correlation function 
and the power spectrum in section [|, demonstrating (by ex- 
ample) that it is possible to set up power spectra yielding 
negative correlation functions without violating any obser- 
vational constraints. Section ^ considers the effects of nega- 
tive correlation functions on the evolution of clustering and 
shows that the existence of anti-correlated regions yields an 
upper bound on the correlation function at scales shortward 
of the scale at which the linear £ goes negative. Finally, we 
use a synthetic NSR (whose form is derived in an Appendix) 
in section H, to demonstrate the above results. 



2 MATHEMATICAL PRELIMINARIES 

We will consider the evolution of a system starting from 
Gaussian initial conditions, with an initial power spectrum 
Pin(fc). The two-point correlation function, £(a, x), is defined 
as the Fourier transform of the power spectrum : 



£(a,x) 



jP (a, k) e lk x 



(27T) 



(2) 



where a is the expansion factor of the Universe (a oc i 2//3 for 
an = 1 matter-dominated Universe). It is usually more 
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convenient to work with the correlation function averaged 
over a sphere of radius x, given by 



[l+£(a,a:)] =l 3 



(9) 



£,(a,x) = I £(a,y)y 2 dy 



(3) 



The power spectrum and the mean correlation function are 
related by the equations: 

3 f°° dk 

£(a,x) — / — P(a,k) [sin(fcx) — kx cos(fcr)] (4) 



and 
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P(a,k) — — I dxx£(a, x) [sin(kx) — kx cos(kx)] . (5) 
3fc Jo 



where I is some other length scale. In the linear regime, 
( < 1 and I » x. At later stages of the evolution, how- 
ever, as f increases, equation (H) shows that the scale x 
becomes smaller and smaller as compared to the scale I. 
Thus, the evolution of clustering at a scale x is determined 
by the original linear power spectrum at the (larger) scale 
I; this suggests that the true non-linear correlation function 
£(a, x) at some scale x can be expressed in terms of the 
linear correlation function £z,(a, Z), evaluated at a different 
scale I. The general solution is expressible as a non- linear 
scaling relation between £(a,x) and £l(o,,x), with I and x 
related by equation (g) ; this can be written as (Nityananda 



In the linear regime, it can be shown that £ ex a — so that 
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£l(o,,x) oc a 2 £ in (a,, x). 
Note that the average two-point correlation function for a 
power spectrum given by P(k) — 5r>(k — k ), where Sn rep- 
resents the Dirac delta function is 
3 

?(0 = » oi-). [sin(fc o ;) - k l cos(fcj)] , (6) 
2-K z l i k 

which is negative over several ranges of scales, crossing zero 
for the first time when tan (k l) = k l, i.e. I ~ 4.493/fc . 
In fact, negative correlation functions always arise for suffi- 
ciently peaked functional forms of the power spectrum. 

The two-point correl ation function satisfies the equation 
of conservation of pairs (Peebles 198C): 



i)t 



dx 







(7) 



where v(a,x) denotes the mean relative velocity of pairs 
at separation x and epoch a. Equation ([?]) can be rewrit- 
ten in terms of the mean 2-point function, £, by defining a 

new d irppnsinnlpss parr vrWity h(n t) = —i, /a'r- this yiplds 



( Nityai anda fc Padnia.iiabliari 1994[) 



(8) 



where we have introduced the variables D = In (1 + £), 
A — In a and X — In x. 

We will now assume h = h [if (a, i.e. h depends on a 
and x only throug h £. This is a standard assumption in the 
current literature (tflamilton et al. 1991 ; |Nitvananda fc Pad- ' 



Pad 



manabl Lan 1994j; |Mo at al. 199"5^ |Padmanabhan 1996a| 
manabl Lan & Engineer 199£ ) and appears reasonably val- 
idated by numerical simulations. Our analysis and results 
depend on this assumption, which may, of course, be looked 
upon as another way of truncating the BBGKY hierarchy. 
In the present work, this assumption will be treated as a ba- 
sic postulate, validated by simulations; we will not address 
the question of its limits of validity. 

Given the above assumption, h = h [<f (a, x)l , one 
can integrate equation J^) to obta in its general solution 



(Nityananda & Padmanabhan 1994), subject to the condi- 
tion that this reduce to the form £ oc a 2 for ( « 1. The 
characteristics of equation (W satisfy the constraint 



£z,(a,Z) = exp 



dz 



h(z)(l + z) 



(10) 



with I = x{l+£{a,x)] 1/3 . 

Given the form of h(£), one can now relate the non- linear 
correlation function to the linear one; this form is usually 
obtained from simulations. However, it can also be shown 



from general theoretical arguments (Padmanabhan 1996a) 
that £(a, x) can be expressed in terms of £l(o.,1) by the 
relations 

f(a,x) oc ^ L {a,l) 3 (1< f< 200, 1 < £l < 5.85) (11) 
Ci(a,0 3/2 (200 < £,5.85 < ft) , 

if we confine ourselves to models with > everywhere. 
More exact fitting functions, valid over th e entire range of 



£ > 0, can be obtained from simulations (Hamilton et al 



199j ; [Peacock fc Dodds 1994J [Mo at al. 199^ )7 for example, 
the functional form given by Hamilton et al. (1991): 



z + 0.358z 3 + 0.02362 6 
1 + 0.0134z 3 + 0.0020z 9 / 2 



(12) 



where z = ft(a, t) = a 2 fin, and fi n is the initial mean corre- 
lation function. 

Thus, given an initial power spectrum Pi n (k), one can deter- 
mine the average correlation function (and, hence, the power 
spectrum, usingequation (B)) at any epoch, by simply using 
equations (Q), Q and the non- linear scaling relation (|l2"|). 

There is, however, one serious problem with these fit- 
ting functions. The original analysis, leading to equations (8) 
and (9) made no assumptions regarding the sign of £l and 
is valid even if £l < at some scales. The fitting functions, 
however, require that the linear mean correlation function 
be everywhere positive. This is partly due to the fact that 
standard power spectra used in cosmology (such as the CDM 
power spectrum and its variants) do have positive definite 
correlation functions and there appeared to be no compul- 
sion to consider the case of negative £. We feel, however, 
that this situation is unsatisfactory for two reasons : 
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(i) The NSR has proved to be a valuable tool in under- 
standing the physics of gravitational clustering and — in 
particular — the transfer of power between two scales. The 
basic question in the study of transfer of power is how a 
sharply peaked power spectrum evolves due to gravitational 
clustering. To answer this question using the NSR, we need 
to generalise it for negative values of the mean correlation 
function. 

(ii) As we will show in section ^, observations cannot rule 
out the existence of negative £ over certain ranges of scales 
in the universe. Since this feature has very important impli- 
cations, it should be taken seriously, without any theoretical 
prejudice. 

It should be noted that attempts have also been made to 
write NSRs directly f or the function A 2 (fc) = fc 3 P(fc)/(27r 2 ^ 



in freq uency space (Peacock & Dodds 1994; Peacock & 



Dodds 1996), rather than in terms of the mean correlation 
function £, as discussed above. In such Fourier space NSRs, 
one writes 



F[AHk L )} 



with the wavenumbers kfjL and &l related by 
k L = [1 + A 2 NL {k NL )]- 1/3 k N L ■ 



(13) 



(14) 



For example, Peacock & Dodds (1994) give the following 
form for the function F (we have put g(£l) — 1 in their 
expression, for the case of an = 1, matter-dominated Uni- 
verse) : 



F(x) = x 



1 + 0.2fix + (Ax 



.a/3 



1 + {[Ax] a /[ll.GSx 1 / 2 ])^ 



V/3 



(15) 



with A = 0.84 and a = /3 = 2. Since A 2 (fc) > 0, k NL > k Ly 
and power is transferred from small to large wavenumbers, 
i. e. from large to small scales. However, no mathematical 
basis presently exists for scaling relations for quantities other 
than the mean correlation function £; this is discussed in 
more detail in the next section. 



3 NSRS FOR DIFFERENT STATISTICAL 
INDICATORS. 

We discuss, in this section, a few important general issues 
related to the role of £ in the context of non-linear scaling 
relations (NSRs). Throughout this paper we will use the 
term NSR to mean a relationship between the exact value 
of some statistical indicator Q (which could be £, £, A = 
fc 3 P(fc)/27r 2 , etc) to its value Qu„ as given by linear theory, 
in the form Q(x) = F[Q Un (l)] with I = (1 + Q) 1/3 x, where 
x and I are two scales and F is a prescribed function. In the 
case of Fourier space quantities like A(fc) or P(k), we will, 
of course, interpret x and I in term of A; -1 for the two scales. 

There are two broad perspectives one can take regard- 
ing such NSRs. In the first approach, NSRs can be treated as 
a set of relations which are of considerable practical utility 



in studying the evolution of clustering. In this approach, 
they are merely convenient sets of approximate rules by 
which non-linear quantities can be obtained from linear ones, 
thereby facilitating a comparison between theory and obser- 
vations. Taken in this spirit, the key issues in this area are 
only the accuracy of the fitting functions F for the NSR, 
the dependences on various parameters and, finally, given 
an NSR for a particular quantity Q, how best to write down 
(purely as fitting functions) NSRs for other statistical indi- 
cators. There are no fundamental issues; NSRs exist entirely 
by accident, but are of use in studying non- linear structure 
formation. 

The above attitude, however, is tantamount to sweeping 
the entire issue under a carpet, without much investigation. 
In fact, it appears that NSRs are indicative of a key (and not 
completely understood) feature of non-linear gravitational 
clustering and hence need to be investigated thoroughly. In 
support of this view, we would like to stress the following : 

(i) Though the original NSR for £ was obtained as a 
fitting function to simulation data by Hamilton et al. (1991), 
the work by Nityananda and Padmanabhan (1995, hereafter 
NP) clearly spelled out its theoretical origins. In particular, 
this work showed that (a) there exists an exact equation 
satisfied by £ whose integral curves lead to the relation Z 3 = 
x 3 (l + f) (b) the only key assumption which is needed to 
obtain the NSR for £ is that h is a function of £ alone. 

(ii) So far, there has been no evidence for any mathe- 
matical foundation for NSRs for other statistical indicators 
like A 2 or even £. In fact, there is no analogue for local dif- 
ferential equations like equation (25) in NP, for other statis- 
tical indicators like A 2 or the power spectrum. Thus, NSRs 
appear to exist only for the mean correlation function £, 
implying that £ enjoys a special status as an indicator of 
clustering. The NSRs given in the literature for A(fe), for 
example, arose purely as an afterthought, guided by the na- 
ture of the NSR for 

(iii) Given the form of £, it is possible to obtain A 2 
and £ by simple analytic procedures. However, if the NSR 
is true for £, it can be shown that such a similar NSR — 
as defined above — cannot be exactly satisfied for A 2 or £. 
While there have been attempts in the literature to obtain 
fitting functions for A 2 ^ in terms of A 2 ^ (Peacock & Dodds 
1994, 1996), along the lines of the NSR for £, such exercises 
do not have the same level of fundamental validity as the 
NSR for the mean correlation function. 

The above distinction is not usually of great relevance 
for the following reasons : First, if the power spectrum is 
a smooth, mildly varying, power law with adequate asymp- 
totic properties, an NSR for £ will lead to similar approx- 
imate relations for A 2 , £, etc. Which of them is used in a 
specific context — assuming that this is all one uses the 
NSR for — is purely a matter of convenience. Further, since 
any fitting function is only approximate, the fact that some 
NSRs are more approximate than others is often not relevant 
in practical situations. 

The situation, however, is very different when one moves 
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away from simple-minded power spectra and considers a 
more generic situation. The original analysis of NP (which 
provided the theoretical foundation for NSRs) did not make 
any assumptions about the nature of the initial mean corre- 
lation function or the power spectrum. Consider, for exam- 
ple, the case of a power spectrum which is sharply peaked 
around some value, such as a Dirac delta function or a Gaus- 
sian. For such a power spectrum, A 2 , £ and £ have very dif- 
ferent shapes. In particular, £ is quite flat at small scales 
and decreases towards large scales, while A is close to zero 
at small scales. Since the shapes of £ and A are different, 
the evolution predicted by using the NSR in real space and 
Fourier space will be totally different. This can be seen very 
clearly in the case of a power spectrum which is a Dirac delta 
function or for any power spectrum with a cutoff at short 
wavelengths (i.e. P(k) = for k > k c , with some finite fc c ). 
It is precisely in such situations that the mathematical foun- 
dations for the two NSRs come into question and, given the 
fact that the NSR for £ has a theoretical basis while that for 
A 2 is only a fitting function, it is the former which must be 
used when the results from the two differ. Numerical simula- 
tions of peaked power spectra (Bagla & Padmanabhan 1997) 
have shown that, even in such cases, power is indeed trans- 
ferred from larger to smaller scales, as indicated by the NSR 
in real space. In fact, such power transfer to scales smaller 
than the cutoff in the power spectrum would be impossible 
if the Fourier space NSR were correct, simply from the fact 
that the power spectrum at a given value of k is a contin- 
uous function of the scale factor (this result is true for any 
power spectrum with a short wavelength cutoff, not merely a 
peaked power spectrum) . While the power spectra of Bagla 
& Padmanabhan (1997) were exponentially damped beyond 
a small range of scales, and hence did not have a sharp cut- 
off, the fact that these simulations did show the transfer of 
power to small scales, as expected from the a;— space NSR, 
also indicates that the correct NSR to use is the one for real 
space, not Fourier space. 

We note finally, that even if theoretical foundations in- 
deed existed for both forms of the NSR and the two yet 
yielded different results for the same initial power spectrum 
/ correlation function, one should take resort to dynami- 
cal tests such as N-body simulations to ascertain which of 
the two is applicable. In fact, such a situation might well 
prove interesting as it would provide information about the 
dynamical correctness of the two scaling relations. 

The strength of the present work is that its conclusions 
do not depend on the form of the NSR but arise entirely 
from its existence. The synthetic NSR of the last section of 
the paper (and the appendix) is used solely to illustrate the 
results. It should be emphasised that we will show analyt- 
ically that the results are independent of the form of the 
NSR. In the spirit of the above discussion, we will work in 
a>space and not in Fourier space during the analysis. 



4 AN ASIDE : OBSERVATIONAL 
CONSTRAINTS 

As mentioned earlier, negative values of £ are not mathemat- 
ically forbidden and, in fact, arise naturally for sufficiently 
peaked forms of the power spectrum. We, will, in this sec- 
tion, briefly review observational constraints to investigate 
whether present-day observations of the correlation function 
or the power spectrum rule out the existence of negative £ 
in our universe. 

At very large scales, I £ lOOO/i -1 Mpc, the CMBR 
anisotropy measurements provide a constraint on the slope 
and amplitude of the power spectrum; £l is probably posi- 
tive in this range. At small scales, l/i" 1 Mpc I 200h~ 1 
Mpc, the galaxy-galaxy correlation function provides the 
shape and amplitude of the baryonic component of the uni- 
verse. This can be converted into providing some handle 
on the underlying dark matter distribution and the results 
are consistent with a positive mean correlation function. 
At still smaller scales, information on the power spectrum 
arises from the study of abundances of bound structures (like 
quasars, damped Lyman-a systems, etc) and it is very likely 
that the mean correlation function is again positive at these 
scales. There is, however, a range of scales (~ 200/i _1 Mpc 
- lOOOftT 1 Mpc) which are not directly probed by present 
observations and one can easily construct power spectra for 
which the correlation function goes negative in this range. Of 
course, while doing this, one has to also ensure that the ad- 
ditional power does not create any observable consequences 
at other scales. 

We demonstrate, next, (by an explicit example) that 
it is possible to construct power spectra which satisfy the 
above observational constraints and yet yield negative £. As 
mentioned above, the observations allow negative values of 
£ over the range 200ft _1 Mpc < x < lOOO/t -1 Mpc. Consider 
the COBE-normalised CDM power spectrum of Bardeen et 
al. (1986) (eq. G3) 



Pcdm(A;) = Ak 



In [I + 2.34g] 
2.34(7 



(16) 



41 -1/2 



[l + 3.89g + (16.1q) 2 + (5.46q) 3 + (6.71g) 

with q — fc/(f2/i 2 Mpc _1 ) (using 9 = 1 in eq. (G3), corre- 
sponding to photons and three types of relativistic neutri- 
nos). The value of A is set by COBE normalisation to be 
about (24h^ 1 Mpc) 4 . A good fit to the average correlation 
function £cdm for the above spectrum is given by Hamilton 
et al. (1991) 



£cdm(x) = 0.51 



In ( 1 + 



1 ln(l + a;/12) 



x/12 

1 + 0.394a; + 0.00316a- 2 



(17) 



1 + 0.142a- + 0.00129z 2 ' 
where x is in f2 _1 /i -2 Mpc and h is the Hubble parameter. 
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0.1 

k, h Mpc~ 



x, h 1 Mpc 



Figure 1. The figure shows the CDM power spectrum (dotted 
line) and the net power spectrum P(k) = Pcdm(^) + P2(k) (solid 
line). It can be seen that the two are different only over a small 
range in fc-space. 



We will work in an Q — 1 Universe and add to the CDM 
spectrum, a peaked spectrum of the form 



P 2 (k) = A'fcV 
the correlation function corresponding to P2(k) is 
72A'\ 3x 2 - 5A 2 



(x 2 + A 2 ) 5 



(18) 



(19) 



which clearly goes negative for x > (5/3) 2 A. 

The net power spectrum P(k) — PcDM{k)+P2{k) (with 
A' ~ (841ft- 1 Mpc) 4 and A = 250 ftMpc" 1 ) is plotted in 
figure (jy) (solid line), along with the CDM spectrum (dotted 
line). It can be seen that the deviation of the net spectrum 
from the CDM spectrum occurs only over an extremely small 
range of k values and this deviation is itself extremely small. 

Figure (^) shows the average correlation function 
corresponding to P(k). It can be seen that £(x) becomes 
negative at x ~ 390ft -1 Mpc and then returns to positive 
values at x ~ 600ft _1 Mpc. In the region probed by COBE 
(lOOOft -1 Mpc - 3000ft" 1 Mpc), £ is indistinguishable from 
the CDM correlation function; this is due to the extremely 
rapid fall-off in £. 2 (x). 

Thus, even broad band power spectra can yield negative 
correlation functions. The example shown in figure (|l|) is just 
one of many possible ways to construct such power spectra, 
which agree with the current observations. It is hence clearly 
important to understand the effects of such anti-correlated 
regions on the evolution. 

We emphasise that the power spectra we have constructed 
will satisfy all observational constraints currently available 



Figure 2. The average correlation function £(x), for the two 
power spectra of figure (1). The correlation function correspond- 
ing to the total power spectrum, P(k), can be seen to go negative 
at x ~ 390ft — 1 Mpc. The two correlation functions approach each 
other beyond x ~ 800ft — 1 Mpc. 



but will introduce some new features in the evolution of 
the universe. The purpose of doing this is only to stress how 
little we know about our universe; we have hence not worried 
about whether there exist suitable physical processes which 
will lead to such a power spectrum. 



5 CONTINUITY ARGUMENTS AND BOUND 
ON f 

We will now consider the evolution of £ with time, at a given 
scale x, and show that the mere existence of a negative value 
of £l, at some scale I, implies an upper limit to the value of 
I which can be reached at any scale smaller than /; further, 
the maximum value £ max at any scale is a function of the 
scale itself and the location Lo, where £j, first crosses zero. 
This effect also results in the correlation function acquiring 
a x~ 3 dependence at late stages of the evolution. 

We begin by noting that the arguments leading to the 
paradigm for transfer of power (see Section 2) made no 
assumption about the sign of £; these arguments, as well 
as equations (|^) and (px|), thus hold good for the negative 
regime as well. Next, the relation 

l = x[l + Z(a,x)] 1/3 , (20) 

which governs the tranfer of power from one scale to another, 
implies that a given value of £ can arise at some scale x only 
by the transfer of power from a single scale I. For | > 0, 
x < I, i.e. power flows to smaller scales, while, if £ < 0, x > I, 
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i. e. power transfers to larger scales. We stress that equation 
( pp| ) made no assumptions about the sign of £ except, of 
course, the obvious condition, (1 + £) > 0. 

Let Lo be the first scale at which £l crosses zero and 
let xo be the scale of interest, with xo < Lo- Initially, in 
the linear regime, f <C 1 , 2:0 is influenced by scales I such 
that (riio- At later epochs, as £ increases, xq is influenced 
by larger and larger scales, upto the scale Lo. Power from 
the scale (infinitesimally larger than) Lo, however, must be 
transferred to scales larger than Lo, since £l < here; 
clearly, this scale (I — Lo + Al) cannot influence xq. How- 
ever, as has already been emphasised, a given value of f can 
arise at the scale xo only by the transfer of power from a 
single scale; this immediately implies that £(xo) can never 
exceed the value £ , given by 
3 



So 



Lo 
x 



(21) 



since this would require power transfer from the scale Lo+Al 
to xo, which is impossible. Similarly, if £l < over some 
range of I, say, from £1 to I2, then no value of £ which satisfies 
the relation 

3 



s 



V 



xo 



1 



(22) 



can be reached at xo, for I' between li and £2- 
The correlation function, however, is a continuous function; 
clearly, then, if some values of £ cannot be reached, £ must be 
bounded from above by the lowest of these values. Thus, the 
maximum value which £ can take at some scale x (x < Lo) 
is given by 

" L,l3 -1 (23) 

where Lo is the smallest scale at which crosses zero. Thus, 
at late times, £ 3> 1, we must have 

-3 



£ OC X 



(24) 



for x smaller than Lo. This is a generic result, independent 
of the exact form of the initial correlation function, except 
for the requirement that it become negative at some scale. 

Note that £(x) oc x~ 3 over some range of scales implies 
that £(2) = over this range (and vice- versa). For example, 
if £(x) m for x > L, then 

i(x>L) = - 



y S(y)dy + y S{v)dy 



(25) 



r A 



(26) 



where A = J y 2 £,(y)dy is a constant. Thus, S( x ) °c x ~ 3 > 
for x > L, if ~ here. The result that the relation 
£(x) oc x~ z spreads to smaller scales at late times can be 
then seen to imply that £(x) becomes negligible at smaller 
and smaller scales at later and later epochs. 

It may appear strange that the mere existence of nega- 
tive £t has far-reaching effects on the evolution. However, it 



must be emphasised that the positive and negative £ regimes 
are not equivalent, due to the lower bound of ( — 1) on £ in 
the negative regime; no such bound exists for £ > 0. This 
implies an asymmetry in the very structure of the equations, 
if one uses £ as a variable. In fact, this also indicates that 
Log[I + £] is probably the correct variable which should be 
used; in this case, the above asymmetry does not arise, since 
Log[l + S\ takes all values from —00 to +00. 

In general, linear theory ceases to be valid when |£| is 
of the order of unity. However, while using the NSR, one 
maps the linearly extrapolated value £l to the actual value 
of £. This poses no mathematical problems for £ > 0, since 
arbitrarily large positive values of linearly extrapolated £l 
are mathematically allowed and map to still higher values of 
actual £ through the NSR; the situation is quite different in 
the negative regime since the actual £ is constrained to be 
greater than (—1). This asymmetry appears to play a role 
in situations in which negative correlations exist. 

We mentioned in section ^ that power transfer to scales 
smaller than the cutoff in the initial power spectrum is im- 
possible for power spectra with sharp cutoffs, if the Fourier 
space NSR is correct. This can be seen by an argument sim- 
ilar to the one discussed above. Consider an initial power 
spectrum with a short wavelength cutoff at some scale k . In 
Fourier space NSRs (see equations ( |l3| ) to (^)), the power 
P(k) at a wavenumber k at some e poch originates at a fixe d 
scale ki, with k\ and k related by ( [Peacock fc Dodds 1994 ) 



k 3 = kf [f + A 2 (fc)] 



(27) 



with A 2 (fc) = k 3 P(k)/(2n 2 ). Since A 2 (fc) > 0, clearly 
ki < k. If we now consider a wavenumber k beyond the cut- 
off in the initial spectrum, with k > k , there exists a range 
of wavenumbers between k and k , where Pi n (k) = and 
the power at k is then forbidden to take values which could 
originate between k and k. Since the power spectrum is a 
continuous function of the scale factor, if P(k) is not allowed 
to take a set of values, it also cannot take values larger than 
the lowest value of this set. Thus, even in Fourier space, it 
can be seen that an upper bound exists on the power spec- 
trum at large wavenumbers (i.e. small scales) k > k , for 
initial power spectra which cut off at some wavenumber k . 
We note that such power spectra do, in fact, give rise to neg- 
ative correlation functions; one could thus have instead used 
the real space NSR to demonstrate the upper bound on £. 
However, there are also other power spectra such as Gaus- 
sians, which do produce negative £ even without sharp cut- 
offs in the initial spectrum; in such cases, the upper bounds 
are not as obvious from a Fourier space analysis but can be 
clearly seen in real space. Of course, as discussed in section 
^, different results may then be obtained for the late-time 
behaviour of P(k) and £ from the two NSRs, especially in 
the case of peaked initial power spectra. In such situations, 
£ is picked out as the statistical indicator to be used in the 
NSR, by the analysis of NP which provides a theoretical 
basis for the real space scaling relations. 
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6 NSR FOR NEGATIVE £ 

We require an NSR for the negative £ regime to explicitly 
illustrate the results of the preceding section. However, the 
form of the NSR given in equation (113) tacitly assumes that 
£i(a, I) > for all I; as mentioned earlier, this requirement 
has no fundamental significance. The above form, however, 
can be derived from theoretical argu ments related to spher - 
ical collapse, in the regime £ > (Padmanabhan 1996a). 



This suggests that the NSR for the case £ < may also be 
obtained by an analysis of a spherically symmetric model, 
with negative correlations at some scale. This will allow us 
to guess at an ansatz for h(£), based on the relation between 
h and the density contrast S for this model. This form for 

will then be used in equation (jl(]) to obtain an NSR in 
the negative £ regime. 

We will use the ansatz, 



<>(«) 



(-!<?< 0) 



(28) 



for ft(£) (see Appendix for motivation). Substituting for h(£) 
in equation (fid), we obtain, after some algebra 



g (1 + g/4) 

(1+|)V3 



1/3 



(29) 



The inverse mapping, giving £ in terms of £l, can be fitted 
to better than 5% by the relation 



a + 0.88^ + 2.086g 



(-i<€<0) 



(30) 



1.798£| - 2.086^ 

Equations (0) and f jjo] ) implicitly determine £(a,x) in terms 
of the linearly evolved correlation function, £l(o,1) for all 
values, positive and negative, of £l. We mention that the 
above NSR, equation (J3C|), will be used solely for the purpose 
of illustrating the results of Section (4), since these results 
are independent of the exact form of the NSR. The true 
mapping between £ and £l in this regime should, of course, 
be determined from N-body simulations. 

We will use the net power spectrum P(k) of figure (|l|) 
(which is allowed observationally) to demonstrate the effects 
of negative £. The results are shown in figures (j3|) and (Q), 
in which £ is plotted against the scale factor, for two scales, 
xi = 150/i _1 Mpc and X2 = 20/i _1 Mpc. It can be seen 
that, i n both cases, £ does not continue to increase with a 



(or, eqi ivalently, with time), but instead flattens out at late 
times and approaches a value of (Lo/xo) 3 — 1 (dotted lines), 
as predicted by the analysis. Here, Lo is the scale at which 
£l first goes negative and xq is the scale under considera- 
tion. We emphasise that the mapping for negative £ has no 
influence on the results; any mapping for this regime would 
have yielded the same results since these stem entirely from 
the existence of negative values of the correlation function. 

Figure (|s|) shows a log-log plot of £ versus x at four 
differe nt epochs; the solid line has a slope of —3. It can 



be seer that more and more scales acquire the ( a i 
behaviour with time. Scales with x Ss 3ft. - Mpc (i.e. £ <; 



x = x, = 150 h 1 Mpc 
L„ ~390 h 1 Mpc 



log [a] 



Figure 3. The figure shows £ as a function of scale factor, at 
x\ = 150/i _1 Mpc, for the power spectrum P(k) discussed in the 
text. The dotted line corresponds to £ oc a~ 2 . 



2x 10 6 ) have not as yet (upto a = 5x 10 s ) been influenced by 
the region with £x, < and will hence only acquire the x~ 3 
slope at later epochs. It should be emphasised that the above 
discussion regarding the upper bound on £(ie) is applicable 
only for scales x which are smaller than the scale Lo at which 
the initial mean correlation function £t first goes negative 
and not for scales x > Lo . We do not discuss the evolution of 
£(a;) for x > Lo here, since this would depend on the exact 
form of the scaling relation in the negative £ regime and we 
presently only possess a rather synthetic NSR in this region. 

Finally, we note that the "asymptotic" behaviour of £ 
does not arise at any particular value of £ but depends only 
on the scales in question, xq and Lo. For example, in figure 
(^), the asymptotic behaviour arises when £ ~ 16 (i.e. in 
the quasi-linear regime) while, in figure (^), this behaviour 
arises when £ ~ 7 x 10 3 (i.e. in the non-linear regime). 

The mean correlation function £ is known to scale as 



£ oc 



and a (in 3 dimensions) in the linear. 



linear and non-l i near regimes of structu re formation ( Hamil- 
Padmanabhan 1996a). However, the exis- 



ton et al. 1991 



quasi- 



fence of negative correlation functions appears to give rise 
to a fourth regime in gravitational clustering, besides the 
above three regimes. Here, £ asymptotically approaches a 
constant value, at any scale shortward of the scale at which 
(,l crosses zero. Clearly, high dynamic range simulations will 
be needed to test this result; this may well prove difficult, or 
even prohibitive, in 3-dimensions. Such high dynamic ranges 
can, however, be attained in 2-D simulations; we note that 



earlie r 2-D simulations (Bagla, Engineer & Padmanabhan 



1997) have shown that £ does not continue to grow as a in 
the non- linear regime (as expected from stable clustering), 
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2 4 6 -2 

log [a] log [x/(h-' Mpc)] 



Figure 4. The figure shows § as a function of scale factor, at 
X2 = 20h~ 1 Mpc, for the power spectrum P(k) discussed in the 
text. The dotted line corresponds to £ oc a~ 2 . 



but begins to flatten out at late times. It is possible that this 
is indicative of the existence of a fourth regime, beyond the 
stable clustering phase. Of course, the present results 
are quite independent of the stable clustering 
ansatz, since they stem solely from the existence 
of the non-linear scaling relation. 

We note finally that the evolution of non-linear struc- 
ture may also be studied using other models; for example, 
one may view the non-linear correlation function as arising 
from t he superpos i tion of the density pro files of non-linear 
halos flSeljak 2000] ; |Peacock fc Smith 2000| ). However, NSRs 
seem to be reasonably validated by numerical simulations; 
we have hence not considered other models, such as the halo 
model, here, but have limited ourselves to attempting to un- 
derstand the scaling relations themselves. 

In summary, the mere existence of negative correlation 
functions have far-reaching effects on the evolution of clus- 
tering. They result in an upper bound on £ at scales smaller 
than the scale at which £l goes negative and also cause the 
relation f oc x~ 3 at these scales, at late times. Negative 
correlation functions are not forbidden observationally; it is 
hence important to use N-body simulations and/or semi- 
analytic arguments to obtain a form for the NSR in this 
regime. 
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7 APPENDIX 

We will use the evolution of a spherically symmetric, com- 
pensated void to obtain some insight into the form of 
Consider an initial density perturbation of the form 



(0 < R < Ei) 



6(R) 







(Ri < R < R 2 ) 
(R 2 < R) , 



(31) 



such that the effects of the under-dense region within Ri are 
"compensated" by the over-dense region between Ri and R 2 . 
This implies that a shell with radius larger than R 2 does not 
feel any excess gravitational force, which would cause it to 
deviate from the smooth background expansion. Such a con- 
figuration can be seen to give rise to negative values of the 
average correlation function. 

We initially investigate the asymptotic behaviour of the h 
function for such a system, and obtain a fitting function for 
h(S) when 5 < 0. These results will be used to argue for 
a form of /i(£), which can be then used to obtain an NSR 
between £ and £l, in the negative £ regime. 



The general equation for the evolution of a spherical shell 
which encloses a mass M can be written as 

d 2 R _ GM 
dt 2 ~ ~~ R 2 ' 
where R(t) is the proper radius of the shell; the mass M is a 
constant, in the absence of shell crossing. The first integral 
of equation is 



(32) 



1 (dR\ 2 GM 
2\lt) =— + E 



(33) 



The constant E gives the total energy of the shell; clearly, 
E > for a shell with radius, R < R 2 while E = for a shell 
with radius, R > R 2 . We will set up the initial conditions 
by assuming that the shell has radius Ri at the initial time, 
ti and moves with the general background expansion at this 
time, with zero peculi ar velocity. The tota l energy of the 
shell is then given by (Padmanabhan 1996b) 



GM 



[sinh 9-9] 



(2£)3/ 2 

The peculiar velocity, v, defined as v = ax, is 

v = -T- lax] — ax = R — HR . 
dt 1 J 

This implies that 



^-HR- 1 



R 

RH 



For an Si = 1, matter dominated Universe, H 
substituting for H and R in equation (B3) yields 



3t_ 
2R 



(36) 

(37) 
(38) 

(2/3*); 

(39) 



We now use equations (|3fp and 
terms of the parameter 6 as 



h = 1 



3 cosh (61/2) 
2 sinh (6»/2) 



sinh l 



(cosh 9-1) 



i) to rewrite the above in 



(40) 



3/2 



Further, the density contrast 8 is given by 
9GMt 2 



1 + 5 = 



2i? 3 



Again replacing for R and t in terms of 9 gives 
_ 9 (sinh 9 - 9f _ n 



(41) 



(42) 



E 



t 



■m- 



where Ki = (l/2)Hf Rf is the initial kinetic energy, 5i is the 
initial density contrast and Hi is the Hubble parameter at 
the initial time; note that Si < 0, so that E > 0. 
The solution to equation (^), with these initial conditions, 
can be written in the parametric form 

GM 
IF 



2 (cosh 9 - If 

In the limit 9 — > oo, 8 — » —1 and h — > —1/2. We thus have 
the interesting result that h — > —1/2 as S — > —1. 
A good fit to the function h(S), using equations ( ^o| ) and 
@, is given by h(S) = 5/3 - S 2 /6 (see figure (|)). The fit 
satisfies both the linear limit (h(S) ~ 5/3, for \8\ <C 1), as 
well as tends to -1/2 as 8 — > —1. The percentage error on 
the fit is less than 10%. 



The preceding analysis of an individual compensated void 
has yielded a fitting function for h(S); we, however, require h 
as a function of £, to find the NSR for the negative £ regime. 
There is - in general - no simple relation between the den- 
sity contrast 8 of a single lump and the correlation function 
f . (Purists will even consider the question of relating | and 
8 meaningless.) Despite this, relations obtained for single 
spherical lumps have proved to be gen eralizable to a descrip- 
tion of statistical quantities like £ ( Pachnaiiabhanl996£ ; 
Padmanabhan 1997t |Mo & White 199ft |Sheth 1998| ; [Engt 
neer, Kanekar fc Padmanabhan 2000| ). In this spirit, and 



R 



[cosh 9-1] 



(35) 



motivated by the relation between h and 8 being quadratic, 
we will attempt an ansatz in which ft is a quadratic function 
of £. In the linear regime, however, h = (2/3)£; this result 
arise s from straightforward application of perturbation the- 
ory (Padmanabhan 1996c) and is independent of the sign 
of £. Next, £ and 8 are both boun ded from below by -1 
and £ tx 8 in the non- linear regime (Engineer, Kanekar & 
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Figure 6. The figure shows the h function, plotted versus density 
contrast, <5, for a compensated void (solid line). The dotted line 
indicates a fit of the form h(8) = a5+b5 2 , with a = 1/3, b = —1/6. 

Padmanabhan 2000). It thus seems reasonable that £ ~ — 1 
when 5 = — 1; this immediately implies that h ~ —1/2 when 
£ fa — 1. We will hence choose 

h(£) = li+l? , (43) 

similar to the form of h(S) and satisfying the condition 
h — > — 1/2 as £ — > — 1. This will be used in equation ( p"c| ) to 
obtain an NSR for the negative f regime. 



© 0000 RAS, MNRAS 000, 000-000 



